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Results for nucleon matrix elements (arising from moments of structure functions) and form factors from a 
mixture of runs using Wilson, clover and overlap fermions (both quenched and unquenched) are presented and 
compared in an effort to explore the size of the chiral 'regime', lattice spacing errors and quenching artefacts. 
While no run covers this whole range of effects the partial results indicate a picture of small lattice spacing errors, 
small quenching effects and only reaching the chiral regime at rather light quark masses. 



1. INTRODUCTION 

Chiral extrapolations of lattice data to the 
physical pion mass, the continuum or a — > limit 
and removal of the 'quenched approximation' re- 
main major sources of systematic uncertainty in 
the determination of hadron matrix elements (see 
eg P). A study of these quantities is necessary to 
understand how QCD binds quarks and gluons to 
form hadronic states and to give an explanation 
about how particle mass and spin arises. Due 
to the importance of proton/neutron scattering 
and DIS experiments most is known about the 
nucleon. In this talk we shall detail our recent 
progress in studying these problems from the lat- 
tice or numerical perspective in particular with 
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regard to the chiral extrapolation. We have com- 
puted v n = (x n_1 ), for n — 2, 3 and 4 related to 
the three lowest unpolarised moments of the F2 
nucleon structure function, which are given by 
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(q = u,d). In this report we shall restrict atten- 
tion to the n = 2 moment only. We shall also 
consider here form factors arising from lepton- 
nucleon scattering, 
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where = 97^9 is the vector current and the 
momentum transfer q = p' — p. 

Specifically we shall consider a variety of ac- 
tions and quark mass ranges, 

1. Wilson fermions at one fixed lattice spac- 
ing, a^ 1 ~ 2.12GeV in the quenched ap- 
proximation at pseudoscalar masses down 
to ~ 310MeV, in order to try to match to 
chiral perturbation theory (%PT). 

2. 0(a)-improved Wilson fermions ('clover 
fermions') in the quenched approximation 
at three lattice spacings a -1 ~ 2.12 - 
3.85GeV with pseudoscalar masses between 
580MeV and 1200MeV to check finite lat- 
tice artefacts. 

3. Unquenched clover fermions at pseu- 
doscalar masses down to ~ 560MeV in order 
to see if there are any discernable quenching 
effects. 

4. Overlap fermions, in the quenched approxi- 
mation at one lattice size a -1 ~ 2.09GeV 
down to pseudoscalar masses of about 
440MeV. These have a chiral symmetry 
even with finite lattice spacing and hence 
have better chiral properties than either 
Wilson or clover fermions. 

Note that the physical pion mass is about ~ 
140MeV and we use r = 0.5fm = (394.6M6V)" 1 
to set the scale. 

What must one achieve to be able to compare 
numerical results with QCD? One hopes that af- 
ter taking the continuum limit, we can match our 
results to known chiral perturbation theory and 
then take the limit m ps — > m„. Practically we 
might thus expect that for a quantity Q of inter- 
est 



Q = FQ(r Q m ps ) + O((a/r ) n ). 



(4) 



F^(r m ps ) describes the (chiral) physics and the 
discretisation errors are 0(a n ) where n = 1 for 
Wilson fermions and n = 2 for clover and over- 
lap fermions. The 0(a) errors may be split into 
terms that remain in the chiral limit, 0(a/A), 
and 0{am q ) ~ O(ar m ps ). While 0(a) improve- 
ment will remove both these terms, near the pion 



mass this second set of lattice errors will be neg- 
ligible anyway, so we do not have to worry about 
this term. Thus the variant we shall try here for 
clover fermions is 

Q = FQ(r m ps ) + d?(a/r ) 2 + d^ar m 2 ps , (5) 

where we shall take and d^ to be constant. 

Naively one might expect a Taylor series ex- 
pansion for to be sufficient, ie 



F?(x) = F^{Q) + c Q x 2 + 



(6) 



Over the last few years expressions for F® have 
been found 



F™»(x) = r a m N (0) 



c? N x 2 



(7) 



F?(x) = v n (0) (1 + c^x 2 \n(x/r A x ) 2 ) + 



In general leading order xPT introduces non- 
analytic terms to the naive Taylor expansion of 
F®, either odd powers, such as x or a; 3 in F™ N 
or logarithmic as in F£ n . Additional terms may 
appear when using the quenched approximation; 
these are shown in the above equations inside 
square brackets. The chiral scale, A x , is usually 
taken to be ~ lGeV. Expressions are known for 
some of the coefficients in the above expansions 
(using for example the results of |2I31 11516] ) and 
rough numerical estimates give 
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-0.66. 
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For quenched QCD an expression for c\ n in terms 
of D and F can be found in (T] giving an approx- 
imate result of 



-0.28. 



2. LATTICE DETAILS 
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2.1. Run parameters 

Details of the statistics and parameter values 
used in the simulations for the nucleon matrix ele- 
ments are given in tables^ d El an dE (where for 
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p 


ft 


Volume 


Confs. 


6.0 


0.1515 


16 3 x 32 


0(980) 


6.0 


0.1530 


16 3 x 32 


0(1130) 


6.0 


0.1550 


16 3 x 32 


0(1360) 


6.0 


0.1550 


24 3 x 32 


0(220) 


6.0 


0.1558 


24 3 x 32 


0(220) 


6.0 


0.1563 


24 3 x 32 


0(220) 


6.0 


0.1563 


32 3 x 48 


0(250) 


6.0 


0.1566 


32 3 x 48 


0(530) 



Table 1 

Parameters for quenched Wilson fermions with Wilson 
glue. (The smallest pseudoscalar mass is ~ 310MeV.) 



P 


ft 


Volume 


Confs. 


6.0 


0.1320 


16 3 x 32 


0(445) 


6.0 


0.1324 


16 3 x 32 


0(560) 


6.0 


0.1333 


16 3 x 32 


0(560) 


6.0 


0.1338 


16 3 x 32 


0(520) 


6.0 


0.1342 


16 3 x 32 


0(735) 


6.2 


0.1333 


24 3 x 48 


O(300) 


6.2 


0.1339 


24 3 x 48 


O(300) 


6.2 


0.1344 


24 3 x 48 


O(300) 


6.2 


0.1349 


24 3 x 48 


0(470) 


6.4 


0.1338 


32 3 x 48 


0(220) 


6.4 


0.1342 


32 3 x 48 


O(120) 


6.4 


0.1346 


32 3 x 48 


0(220) 


6.4 


0.1350 


32 3 x 48 


0(320) 


6.4 


0.1353 


32 3 x 64 


0(260) 



Table 2 

Parameters for quenched clover fermions with Wilson glue. 

Wilson or clover fermions, am q = (1/k— 1/k c )/2 
with k c to be determined). Note that these refer 
to the number of configurations used for finding 
the matrix elements; for the nucleon/pion mass 
sometimes a higher statistic was used. 

2.2. Hadron masses 

To give some idea of where these parameter 
values lie, in Fig. ^ we plot the various nucleon 
masses against the square of the pion mass. Also 
shown for comparison is the position of a hypo- 
thetical ss pseudoscalar meson and the physical 
pion, as well as the chiral limit. The data lie in 
a narrow band tending to the physical nucleon 
(denoted by a star in the plot). In general there 
do not appear to be large discretisation effects. 
For the quenched data (in particular) taking the 



p 


^sea 


Volume 


Trajs. 


Group 


5.2 
5.2 

5.2 


0.1342 
0.1350 
0.1355 


16 3 x 32 
16 3 x 32 
16 3 x 32 


5000 
8000 
8000 


QCDSF 

UKQCD 

UKQCD 


5.25 
5.25 
5.25 


0.1346 
0.1352 
0.13575 


16 3 x 32 
16 3 x 32 
24 3 x 48 


2000 
8000 
1000 


QCDSF 

UKQCD 

QCDSF 


5.26 


0.1345 


16 3 x 32 


4000 


UKQCD 


5.29 
5.29 
5.29 


0.1340 
0.1350 
0.1355 


16 3 x 32 
16 3 x 32 
24 3 x 48 


4000 
5000 
2000 


UKQCD 

QCDSF 

QCDSF 


5.4 
Table 2 


0.1350 


24 3 x 48 


1000 


QCDSF 



Parameters for unquenched clover fermions with Wilson 
glue. 



p 


am q 


Volume 


Confs. 


8.45 


0.140 


16 3 x 32 


O(50) 


8.45 


0.098 


16 3 x 32 


O(50) 


8.45 


0.056 


16 3 x 32 


O(50) 


8.45 


0.028 


16 3 x 32 


O(50) 



Table 4 

Parameters for quenched overlap fermions with (tadpole 
improved) Liischer-Weisz glue. 

y-axis to be oc m 2 N seems phenomenologically to 
lead to linear lines (although possibly the lightest 
point in a data set might be bending slightly up- 
wards). A pseudoscalar mass of ~ 400MeV cor- 
responds to (romps) 2 ~ 1.0 while m ps ~ 600MeV 
gives (rom ps ) 2 ~ 2.3. Most of the lattice data lies 
above the region m ps ~ 600MeV, with the ex- 
ception of the quenched Wilson fermion results. 
In what region might the leading order xPT be 
applicable? In Fig. ^ we a ls° plot a representa- 
tive curve (for the quenched nucleon mass). A 
breakdown is seen at about (rom ps ) 2 ~ 1 or 
m ps ~ 400MeV, as the gradient of the xPT result 
is too steep. Higher order xPT might improve the 
situation |8I9| . 

2.3. Nucleon matrix elements 

There are two distinct steps in determining nu- 
cleon matrix elements. First the bare matrix ele- 
ment must be determined from the ratio of three 
point nucleon-operator-nucleon correlation func- 
tions to two point nuclcon-nucleon correlation 
functions. Secondly this matrix element must 
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■ P=6.00 24x32 Quen, Wilson 
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Fi gure 1. (ro m Jv) 2 versus (romps) 2 for the four sets under consideration. The experimental nucleon mass is denoted by a 



star. 



The chiral limit m,p S = is shown as a short dashed line, while the physical pion mass is denoted by the long dashed 



line. Also shown as a dot-dashed line is the mass of a hypothetical ss meson calculated as ■ 
line uses the numerical values from eq. J5t . 



V2m K . The (quenched) xPT 



be renormalised. Techniques for both steps are 
standard, eg |10| . In general, one calculates non- 
singlet, NS, matrix elements because then the 
difficult to compute one-quark-line-disconnected 
part of the matrix element cancels (although for 
the vector current this additional piece is proba- 
bly very small, so we shall ignore this point 
here). 

For all the data sets for V2 and the vector 
current the appropriate renormalisation constant 
is known to one loop perturbation theory (for 
overlap fermions the results are given in |12I13| ). 
Either one can tadpole improve, TI, the re- 
sult, for example TRB-PT, |14I15| or use a non- 
perturbative method, NP, (For our variant 

of the method see HZ|.) For overlap fermions 
and some of the unqucnched clover results this 
non-perturbative method has not yet been imple- 
mented, so we must use TFPT. As we wish to 
compare all our results, for consistency we shall 
present them all here using TI-PT. The fact that 



the vector current is conserved however allows 
a more direct NP determination of Zy, |11I12| . 
For overlap fermions in particular this appears to 
indicate that there are significant differences be- 
tween perturbative based results and NP results; 
however for Wilson or clover results there is of- 
ten only a couple of percent (or less) difference 
between TI and NP results. 

3. MATRIX ELEMENT RESULTS 
3.1. uf 3 

We now show a selection of results. From the 
quenched clover fermion data we shall first check 
for lattice spacing effects. In Fig. [21 we show the 
chiral and continuum extrapolation for v-ib (the 
'6' just denotes the lattice represention used) us- 
ing the formula given in eq. JHJl, together with a 
linear function in the quark mass, eq. ©• The 
upper graph represents the chiral physics and 
the lower graph represents the discretisation er- 
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Figure 2. u™ S; 

of 2GeV versus (romps) 2 , upper picture, and vi, 
(F^ 2 (0) + c(r m pa ) 2 ) — d m aT a m? vs against (a/r ) 2 , lower 
picture, for quenched clover fermions. Filled circles, 
squares and diamonds represent /3 = 6.0, 6.2 and 6.4 re- 
spectively. The MRS phenomenological value of ?; 2 ' 
is represented by a star. 



rors. First we note that there seem to be very 
small a discretisation errors and indeed the nu- 
merical effect of additional operators needed to 
ensure O(a) improvement seems to be negligi- 
ble, |15| . Secondly a linear chiral extrapolation 
seems to describe the data adequately. (Again, 
as for the nucleon mass, we might expect 
to be only valid for m ps < 400MeV.) A result of 
v^(2GeV) = 0.25(1) is found. This is to be com- 
pared to the phenomenological result of ~ 0.18, 
so we see about a 40% discrepency. 

Turning now to the unquenched data, in Fig. 13 
we show the results. Again similar conclusions 
hold as for the quenched case, although it should 
be noted that we are not so close to the continuum 
limit. Wc find vf^(2GeV) = 0.27(2). So at least 
in this quark mass range quenching effects seem 
to be small. 

Is the use of Wilson or clover fermions bad for 
the investigation of chiral properties and the chi- 
ral limit? Overlap fermions having a chiral in- 
variance on the lattice in the massless limit, are 
in principle better. In Fig.0|we show preliminary 
results for vi using quenched overlap fermions, 
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Figure 3. Same as Fig.[5]for unquenched clover fermions. 
Filled circles, squares, up triangle, diamonds and down 
triangle represent = 5.20, 5.25, 5.26, 5.29 and 5.40 re- 
spectively. 

|12| . Again we see the linear behaviour in the 
quark mass confirmed. (The result is in surpris- 
ingly good agreement with the phenomenological 
value. At present we have to use the TI per- 
turbation result for the renormalisation constant. 
As noted previously when comparing a NP eval- 
uation of Zy with a TI perturbative evaluation 
the differences are greater than in the Wilson or 
clover case. So perhaps a NP evaluation of the 
renormalisation constant will give it a large value 
and so lead to a larger result for V2, more in line 
with other values presented here.) 

Finally we investigate the chiral limit, using 
light Wilson quarks on large lattices. In Fig. [3] we 
show the results for (roTO ps ) 2 between about 0.5 
and 2.5 or m ps - 300 - 600MeV. First it seems 
that the numerical values are about the same as 
for the clover quenched case, which might indicate 
that finite lattice effects are again small. All the 
results are again rather flat (and this continues 
to heavy quark masses |10I14| ). Also shown is a 
linear fit from eq. © together with xPT results, 
eq. @ , starting from the MRS phenomenological 
value and using the quenched result in eq. (fTUf> 
for chiral scale lGeV. (Using larger values of c\ 2 
leads to a steeper slope, while using a smaller 
value for A x reduces the slope.) For the present 



6 



0.4 



0.3 



0.2 



0.1 




0.0 



2.0 



4.0 6.0 

< r o m ps ) 2 



10.0 



Figure 4. u" s ' MS (2GeV) using quenched overlap 
fcrmions at fi = 8.45. 

numerical results to reach the phenomenological 
value seems difficult, as down to ~ 300MeV they 
seem to be rather linear. Of course we would 
not necessarily expect complete agreement with 
the phenomenological result, however past expe- 
rience with the quenched approximation does sug- 
gest that it can lead to values only a few percent 
away from experimental values. 

3.2. Form Factors 

Finally we present our results for electromag- 
netic form factors, as given in eq. 0- We have 
i*i(0) = 1 (charge conservation leading to a NP 
evaluation of Zy on the lattice) andF 2 (0) = fj,— 1, 
the anomalous magnetic moment in magnetons. 
It is usual to give the results for the Sachs form 
factors 



G e (q 2 ) 
G m (q 2 ) 



Fi(q 2 ) + 



{2m N y 
= F^q 2 ) + F 2 (q 2 ) . 



F2(q 2 ). 



(11) 



Previous experimental results give phenomeno- 
logical dipole fits, for the proton, p or neutron, 

71, 



Gl(q 2 ) 



G p m (l 2 ) 



G? H (q 2 ) 



(12) 



0.35 



0.30 



0.25 



0.20 - 



0.15 



Figure 5. v^' MS (2GeV) using quenched Wilson 
fermions at fi = 6.0. The 16 3 X 32 results are shown as 
circles, the 24 3 X 32 results as squares and the 32 3 X 48 
results as diamonds. A linear fit, eq. J^J, is also shown as a 
dashed line. The (quenched) xPT une uses the numerical 
values from eq. 1101 . 



and 




G n e (q 2 ) 



0. 



(13) 



with m v ~ 0.82GeV, fi p ~ 2.79 and fi n 1.91. 

(Present experimental results indicate however 
that, JH|, GP(q 2 )/G P n (q 2 ) is decreasing with in- 
creasing — q 2 .) In Fig. we show G v e and G p m for 
the proton and in Fig. f° r the neutron. Note 
that we have made simple linear chiral extrapo- 
lations here, see eq. (0, but at present we can 
do little more than this, as in particular on the 
large 32 3 x 48 lattice, we only have two quark 
mass values. The technique used here is the same 
as described in ^) 20 . Previous lattice calcula- 
tions described in these papers seem to show little 
or no lattice artifacts. Due to the larger lattice 
size, we are now able to go to smaller momentum 
transfer -q 2 mm - (27r/32a) 2 - 0.17GeV 2 than 
was previously possible. We find that the lattice 
result tracks the experimental values reasonably 
well, although with rather large errors. Finally 
we note that as usual the electric form factor of 
the neutron remains difficult to measure on the 
lattice (see for a recent attempt). 
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Fi gure 6. Ge and G^ for the proton as functions of — q 
using quenched Wilson fermions at /3 = 6.0. The previous 
experimental results are given by the open circles. 
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